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The anomalous properties of water in the supercooled state are numerous and well-known. Par- 
ticularly striking are the strong changes in dynamic properties that appear to display divergences 
at temperatures close to - but beyond - the lowest temperatures attainable either experimentally 
or in computer simulations. Recent work on slow or glassy dynamics in water suggests analogies 
with simple liquids not previously appreciated, and at the same time highlights some aspects that 
remain peculiar to water. A comparison of the behavior of water with normal liquids, with respect 
to its dynamic, thermodynamic and structural changes in the supercooled regime is made by ana- 
lyzing, via computer simulations, the properties of local potential energy minima sampled by water 
in supercooled temperatures and pressures. 



I. INTRODUCTION 

Water is among the most abundant substances on 
earth, and also among the most familiar, as it forms a 
significant part of the natural environment, and of living 
organisms. In addition to these obvious reasons, the con- 
tinued study of water by physical scientists stems from 
the many remarkable and peculiar properties it possesses, 
some of which are part of the popular lore (such as the 
fact that liquid water is heavier than ice near ambient 
conditions). In the crystalline state, water exists in as 
many as twelve, possibly fourteen, forms of ice. Echoes 
of this polymorphism are to be found in the amorphous 
state as well. When prepared as an amorphous solid, or 
glass, by ultrafast cooling or by depositing vapor on a 
cold substrate, water is found in two forms, a high den- 
sity and a low density form |l[ . 

A substantial fraction of the study of liquid water has 
focussed on its properties in the supercooled state [§, 
where its unusual properties are amplified. Analyzing a 
range of dynamic and thermodynamic quantities. Speedy 
and Angell showed that these quantities appeared to 
diverge at a finite temperature in a power law fashion, 
leading to the hypothesis of a thermodynamic singularity 
at Tg = 228K. Computer simulation studies and anal- 
ysis of simplified models have in recent years led to the 
proposal of two other possibilities, namely the existence 
of a liquid- liquid critical point at low temperatures [Q, 
and the scenario wherein no thermodynamic singulari- 
ties need be invoked to explain the anomalies observed 
in water 

However, there has been a renewed focus recently on 
the idea that regardless of a thermodynamic singularity, 
strong changes in both the thermodynamic and dynamic 



properties of supercooled water may be expected in the 
vicinity of Evidence from computer simula- 

tions suggests that the strong temperature dependence 
of dynamical quantities, and the power law tempera- 
ture dependence, may be explained as manifestations 
of an avoided dynamical singularity described by the 
mode coupling theory of the glass transition Anal- 
ysis of experimental data appears to indicate that 
the rapid changes in the entropy of supercooled water 
must cross over to slower changes in the vicinity of T,, 
and correspondingly, the heat capacity must display a 
maximum. It has been argued that upon cross- 

ing this temperature, water changes character from be- 
ing a very "fragile" liquid at higher temperatures, to a 
"strong" liquid at lower temperatures Specifically, 
the stronger than Arrhenius temperature dependence of 
relaxation times seen at temperatures higher than Tg is 
argued to cross over to Arrhenius dependence at lower 
temperatures (the experimental situation is ambiguous 
at present, with measurements in Ref. [ pl| arguing in 
favor of continuous change down to the glass transition 
temperature ^ 135K). 

Among the crucial issues one must address in under- 
standing the above possibilities is the nature of structural 
change that takes place as liquid water is supercooled. In 
particular, the nature of structural change that accompa- 
nies a fragile to strong crossover must be understood. In 
this paper we present a preliminary report of investiga- 
tions to address these questions, which focus on analyzing 
the properties of local potential energy minima or "inher- 
ent structures" sampled by a simulation model of water. 
Similar investigations have recently proved a fruitful ap- 
proach to studying changes in liquid state properties un- 
der supercooling of model atomic liquids as well as liquid 
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water |12|-|17|. In particular, we focus on the comparison 
of changes in dynamics with changes in the energies of 
the inherent structures sampled, and structural change as 
evidenced by changes in the nearest neighbor geometri- 
cal properties. A specific structural feature we consider is 
the fraction of water molecules that are four-coordinated 
and those with five or greater other molecules in their first 
neighbor shell. These latter "bifurcated bond" arrange- 
ments have previously been shown to facilitate structural 
rearrangement and hence faster dynamics [ p^ . 
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FIG. 1. Relaxation times: (a) Shown on a logarithmic 
scale vs. inverse temperature (Arrhenius plot). Deviations 
from Arrhenius behavior is observed at low temperatures, (b) 
Transformed to yield a constant value if Arrhenius behavior 
is obeyed. 



II. SIMULATION DETAILS 

We perform molecular dynamics (MD) simulations of 
216 water molecules interacting via the SPC/E pair po- 
tential The simulations we describe here are for a 
fixed density of p ~ l.Og/cm^, for temperatures ranging 
from 210/^ to 700K. For T < 300 K, we simulate two 
independent systems to improve statistics, as the long re- 
laxation times make time averaging more difficult. Full 
details of the simulation protocol used are described in 
p9[ . At the studied density, the mode coupling singular 
temperature, Tc has been estimated to be Tc = 19AK± 1 

Inherent structures are obtained, with a sample of 100 
equilibrated liquid configurations as starting configura- 
tions, by performing a local minimization of the potential 
energy, using conjugate gradient minimization. The iter- 
ative procedure for minimizing the energy is performed 
until changes in potential energy per iteration are no 



more than lQ~^^kJ /mol. A normal mode analysis con- 
firms that the potential energy changes with positive cur- 
vature along all directions away from the minimum en- 
ergy configurations so obtained. 

We calculate the oxygen-oxygen pair correlation func- 
tion for both the equilibrium liquid configurations, and 
the inherent structures. The integrated intensity under 
the first peak of the pair correlation function (upto the 
cutoff radius O.Slnm) yields the coordination number. 
Dynamics is probed by calculating the self intermediate 
scattering function at wavevector 18.55 nm~^, i.e. the 
Fourier transform of the self van Hove self correlation 
function 
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Stretched exponential fits of the self intermediate scat- 
tering function yield relaxation times r as a function of 
temperature. Average values of the inherent structure 
potential energy and pressure are also calculated. 
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FIG. 2. (a) Average inherent structure energies, (b) Equi- 
librium averaged potential energies. 



III. RESULTS 

Figure 1(a) shows the relaxation times plotted on a 
logarithmic scale against inverse temperature. Clear de- 
viations from the high temperature Arrhenius behavior. 



T = Taexp{~E/kBT) 
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(where E is a. constant) is observed for T < 2757^. This 
is more clearly seen by plotting T ln{T/To) (where tq is 
obtained from an Arrhenius fit of the highest 5 temper- 
atures), as shown in Fig. 1(b), which shows strong devi- 
ation of T ln{T/To) below T ~ 275K indicating that the 
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temperature dependence of relaxation times at low tem- 
peratures can no longer be described by the Arrhenius 
form. 

Figure 2 (a) shows the average inherent structure en- 
ergies as a function of temperature. For comparison, the 
equilibrium potential energy is shown in Figure 2(b). In- 
herent structure energies change continuously, with mod- 
est temperature dependence at high temperatures to sub- 
stantial temperature dependence at lower temperatures. 
In particular, in the range of temperatures where the re- 
laxation times begin to display non-Arrhenius behavior, 
the inherent structure energies show considerable tem- 
perature dependence. This behavior is analogous to that 
found for a model atomic fluid | |T^ . 
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FIG. 3. Oxygen-oxygen pair correlation function, for (a) 
equilibrated liquid configurations, and (b) inherent structures. 
Upon lowering temperature, the first peak in both cases be- 
comes sharper, and the intensity between the first and sec- 
ond peaks decreases. The smaller changes seen in the case 
of the inherent structures oflFers an estimate of that part of 
the change due to configurational change upon cooling, as 
opposed to thermal effects. 

Figure 3 shows the oxygen-oxygen pair correlation 
function, (a) for the equilibrated liquid and (b) for the in- 
herent structures. In contrast to the model atomic liquid 
studied in ||l^ , |l3|| , both pair correlation functions show 
marked temperature dependence, at fixed density. In 
both Fig.s 3(a) and 3(b), the first peak of the pair cor- 
relation function becomes sharper upon decreasing the 
temperature. Further, there is a systematic reduction of 
the intensity between the first and second neighbor peaks. 
The integrated intensity under the first peak, which gives 
the average number of neighbor molecules in the first co- 
ordinate shell, approaches the value of 4 from higher val- 
ues, as the temperature decreases. This implies that the 
configurations sampled by the liquid come closer to that 
of a four-coordinated random network. 
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FIG. 4. Distribution of the coordination number of 
molecules: (a) Equilibrated configurations, and (b) Inherent 
structures. In both cases, the distribution becomes more nar- 
rowly distributed around four-coordination. 

To quantify this further. Fig. 4 shows the histogram 
of the fraction of molecules with a given coordination 
number. For equilibrated configurations, the histogram 
changes from a rather broad one to one that is peaked 
around the value 4 as the temperature decreases. The 
same trend is visible for the inherent structures, although 
even at high temperatures, the distribution is quite nar- 
rowly peaked around the value 4. Such a comparison 
permits us to make a separation between deviations from 
four-coordination arising from thermal agitation pO| , and 
that arising from configurational change. 
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FIG. 5. Fraction of four coordinated water molecules: (a) 
Equlibrated configurations, and (b) Inherent structures 

From data in Fig. 3 and 4, we calculate the fraction 
of molecules that are four-coordinated, and those that 
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have five or higher coordination, i. e. which have bifur- 
cated bonds. Fig. 5 shows the temperature dependence 
of the fraction of four-coordinated water molecules. For 
both the equilibrated configurations and inherent struc- 
tures, this fraction approaches one as T — > Tc- As with 
deviations from Arrhenius behavior, the range of temper- 
atures displaying substantial increase in this fraction also 
shows substantial temperature dependence of the average 
inherent structure energies. 

Figure 6 shows the fraction of molecules with five or 
higher coordination. Complementary to the variation of 
the fraction of four coordinated molecules, this fraction 
approaches values very close to zero as Tc is approached. 



the rate of change in entropy of the liquid would change 
substantially near the crossover temperature. Similarly, 
because of the significant temperature dependence of the 
fraction of bifurcated bonds - which facilitate structural 
rearrangement - above the crossover temperature, and 
the relative constancy below, the temperature depen- 
dence of the dynamical properties may also be expected 
to show a corresponding crossover. Further work is in 
progress to strenghten these notions. 

We acknowledge useful discussions with Prof. C. A. 
Angell. FS acknowledges INFM-PRA and MURST- 
PRIN, and HES acknowledges NSF, for financial support. 
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Fraction of five or higher coordinated water 
(a) Equlibrated configurations, and (b) Inherent 



IV. SUMMARY 

We have presented simulation results that demonstrate 
the significant correlations between relaxation times, av- 
erage energies of local potential energy minima sampled, 
and structural features for a simulation model of water. 
In particular, as the mode coupling temperature Tc es- 
timated for this model liquid at the studied density is 
approached, the structure of the liquid appears to ap- 
proach that of a four-coordinated network, free of defects 
to the extent permissible by the bulk density of the liq- 
uid. Such a structural change can potentially explain the 
speculated dramatic changes both in the dynamics and 
thermodynamic properties of water across this crossover 
temperature. Below the crossover temperature, no fur- 
ther structural arrangement may be expected, and the 
configurational entropy of the liquid may become 'frozen 
in' at the value that corresponds to the random tetra- 
hedral network (plus the residual defects that may be 
present at concentrations varying with density). Thus 



[1 
[2: 

[s: 
[4: 
[s: 
[6: 

[7 
[8 

[9: 
[10: 

[11 

[12 
[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 



Email: sastry@jncasr.ac.in 

O. Mishima, L. D. Calvert and E. Whalley, Nature 288, 
76 (1985). 

C. A. Angell, in Water: A Comprehensive Treatise, F. 
Franks, ed. (Plenum Press, New York, 1972) vol. 7. 
R. J. Speedy and C. A. AngeU, J. Chem. Phys. 65, 851 
(1976). 

P. H. Poole, F. Sciortino, U. Essmann and H. E. Stanley, 
Nature 360, 324 (1992). 

S. Sastry, P. G. Debenedetti, F. Sciortino and H. E. Stan- 
ley, Phys. Rev. £53, 6144 (1996). 

F. Sciortino, P. Gallo, P. Tartaglia and S.-H. Chen, Phys. 
Rev. E 54, 6331 (1996). 

K. Ito, C. T. Moynihan and C. A. AngeU, Nature 398, 

492 (1999); S. Sastry, Nature 398, 467 (1999). 

F. Starr, C. A. AngeU, R. J. Speedy and H. E. Stanley, 

(preprint). 

C. A. Angell, J. Non-Cryst. Solids 131-133, 13 (1991). 
F. Sciortino, A. Geiger, and H. E. Stanley, Nature 354, 
218 (1991). 

R. S. Smith and B. D. Kay, Nature 398, 788 (1999); K.P. 
Stevenson, G. A. Kimmel, Z. Dohnalek, R. S. Smith, and 
B. D. Kay, Science 283,1505 (1999). 
S. Sastry, P. G. Debenedetti and F.H. Stillinger, Nature 
393, 554 (1998). 

F. Sciortino, W. Kob and P. Tartaglia, Phys. Rev. Lett. 
83, 3214 (1997). 

A. Scala, F. W. Starr, E. La Nave, F. Sciortino and H. 



E. Stanley, Eond-mat/9908301 



L. Angelani, G. Parisi, G. Ruocco and G. Viliani, Phys 
Rev. Lett. 81 4648 (1998). 

S. Buechner and A. Heuer, :ond-mat /990627^ ; cond- 



mat/9906280 



C. J. Roberts, P. G. Debenedetti and F. H. Stilhnger 
(preprint). 

H. J. C. Berendsen, J. R. Grigera, and T. P. Stroatsma, 
J. Phys. Chem. 91, 6269 (1987). 

F.W. Starr, F. Sciortino, and H.E. Stanley, Phys. Rev. E 
60, 6757 (1999). 

T. A. Weber and F. H. Stillinger, J. Chem. Phys. 87, 
3252 (1987). 



4 



